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The temperature dependence below Tc of the Uneshape of optical phonons of dif- 
ferent symmetry as seen in Raman scattering is investigated for superconductors with 
anisotropic energy gaps. It is shown that the symmetry of the electron-phonon vertex 
produces non-trivial couplings to an anisotropic energy gap which leads to unique 
changes in the phonon lineshape for phonons of different symmetry. The phonon 
lineshape is calculated in detail for Big and Aig phonons in a superconductor with 
d^2_y2 pairing symmetry. The role of satellites peaks generated by the electron- 
phonon coupling are also addressed. The theory accounts for the substantial phonon 
narrowing of the Big phonon, while narrowing of the Aig phonon which is indistin- 
guishable from the normal state is shown, in agreement with recent measurements 
on Bi2Sr2CaCu208 and YBa2Cu307. 
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I. INTRODUCTION 



Optical phonons observed via Raman scattering have provided a large amount of in- 
formation concerning the energy gap in high-Tc superconductors [|l|, and there have been 
attempts to describe the changes in the phonon lineshapes below Tc in s— wave and 
(i— wave superconductors 0. It is believed that the changes in the phonon lineshape below 
Tc are due in part to changes in the phonon self-energy resulting from coupling between 
phonons and quasiparticles. It has been argued that if the optical phonon has a frequency 
below the pair threshold energy 2A, then the phonon's linewidth decreases (narrows) and its 
frequency renormalizes to lower frequencies (softens) as the quasiparticles become frozen out. 
However, for a phonon near 2A, the linewidth is predicted to grow due to the enhancement 
of the density of states at the gap edge and there can be either pronounced phonon softening 
or hardening depending on which side of the threshold the phonon is located. This simple 
picture has been employed to determine the position of 2A in the cuprate superconductors. 

However, this simple analysis applied to the cuprate systems has revealed that the above 
picture is a bit misleading. The above scenario has yielded a value for the energy gap that 
is different for different types of optical phonons and is thus symmetry- dependent. For the 
case of the Bi 2:2:1:2 system, where very clean surfaces can be obtained, a low frequency 
phonon which transforms according to Aig symmetry (located at 464 cm~^, connected with 
the bridging Oxygen vibrations) shows a downward frequency shift (softening) while no sub- 
stantial linewidth change from the normal state can be resolved from the data However, 
the Big phonon (285 cm~^, connected with the antisymmetric out of plane 0(2) and 0(3) 
vibrations in the Cu-0 plane) on the contrary shows a small frequency softening but a sub- 
stantial linewidth narrowing below T^ 0. Similar behavior is seen in the YBCO systems 
0, where such a large difference in behavior between the Aig and Big phonons in part led 
the authors of Ref. ^ to suggest that these two phonons interact with different electronic 
systems. 

There has been no satisfactory theoretical explanation for the behavior of the different 
phonons. The main problem in addressing these experiments with the existing theories 
concerns the lack of attention paid to the symmetry dependence of the optical phonons. 
However, this symmetry dependence can be an important tool to uniquely determine the 
k— dependence of the energy gap around the Fermi surface. It has been shown that the 
electronic contribution to Raman scattering can provide a large amount of polarization 
(symmetry) dependent information that allows for a stringent test to made to determine the 
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actual symmetry of the energy gap in superconductors 0. It was shown that the couphng 
between the Raman vertex and an anisotropic gap leads to symmetry dependent spectra, 
with peak positions and low frequency and temperature behavior dependent on polarization 
orientations. These changes in the spectra allow for a direct determination of | A(k) |. 
Good agreement with the electronic Raman spectra taken on very clean BSCCO surfaces 
was obtained using a gap which was predominantly or entirely of d^2_y2 symmetry, where 
the peak position and the low frequency behavior of the spectra could be straightforwardly 
accounted for. The symmetry dependence of the data led to the conclusion that the gap 
must be predominantly of Big character. Since the phonon self-energy is very similar to the 
electronic Raman density response, the same type of analysis for the electronic contribution 
to Raman scattering can be made to the phonons as well, leading to a further check on the 
predictions recently made concerning the energy gap in the cuprate materials. 

We propose an alternative explanation for the symmetry dependence of the Raman shifts 
based upon nontrivial couplings of phonons of different symmetry with an anisotropic energy 
gap. Close attention will be paid to the role of the electron-phonon vertex, and consequences 
of its k-dependence will be addressed. Most importantly, it is shown that the lineshape is 
polarization dependent for anisotropic superconductors and different dependences on tem- 
perature can be used to determine not only the magnitude but the symmetry dependence 
of the energy gap. Moreover, it is shown that the peak of the self energy can be located 
at frequencies below 2Amax for certain polarizations which have an symmetry orthogonal 
to that of the energy gap. Thus if the symmetry of the phonon is neglected, values of the 
energy gap inferred from changes in the phonon lineshape using an isotropic s-wave theory 
will be underestimated. In particular, the phonon spectral function for a superconductor 
with d^2_y2 symmetry is examined and a comparison is made with experimental data on 
both the Big and the Aig phonons in BSCCO and YBCO. It is shown that satisfactory 
agreement can be obtained which reconciles the differences between the Aig and Big phonon 
lineshapes. 



II. PHONON SPECTRAL FUNCTION 



The phonon spectral function is given by 
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where uq is the optical phonon frequency and S', S" are the real and imaginary parts of the 
phonon self-energy, respectively. The real part of the self energy renormalizes the position 
of the phonon, while the imaginary part governs the linewidth. The interaction of optical 
phonons and electrons can be simply written as 

He-ph= ^k(q)cLq,<xCk,a(fei,7 + feq,7), (2) 

k,q,7,(7 

where g^{ci) is the matrix element for scattering an electron from k ^ k — q, and 6q ,^, 6]^^ 
are the field operators for phonons of branch 7. The details of the scattering matrix elements 
depend on the nature of the mechanism of the electron-phonon coupling and the symmetry of 
the lattice vibration. In this paper we only consider the symmetry of the matrix element and 
leave a treament of the mechanism and magnitude of the coupling for future consideration 

§• 

We take the k dependence along the Fermi surface of the vertex into account by expanding 
in terms of Fermi surface harmonics $ for small q, 

9l = T.9im^), (3) 

L 

where the index L indicates the order of polynomial that transforms according to the — th 
representation of the point group of the crystal. For cylindrical Fermi surfaces, L can 
be replaced by azimuthal quantum numbers. The symmetry of the optical phonon enters 
into the matrix elements g^. The matrix elements for the phonons accessible to in-plane 
polarizations are given for a cylindrical (2-D) Fermi surface in terms of azimuthal angle (p 
as 

^7f''=^?f=2V2cos(20) + ... (4) 

where we have dropped higher order terms, arguing that they are more anisotropic than 
the terms considering here and will hence be of minor importance. The L = 2 term for 
the Aig channel which is present for z dispersion is absent here and the L = 4 term is 
the first anisotropic term in the series in this case Also, since there is no dispersion in 
the z direction in this case, there are no contributions to the Eg channels. Consequences 
of the Fermi surface and the resulting response functions are considered in a forthcoming 
publication 0, and thus for our purposes we will confine our attention to only cylindrical 
Fermi surfaces. 
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The form of the e — ph interaction Eq. (2) is similar to the electronic contribution to 
Raman scattering in the case of non-resonant scattering with the replacement of the effective 
Raman vertex by the e — ph coupling vertex. Thus we can proceed along the lines recently 
taken for the case of the electronic Raman scattering , where it was shown that the Raman 
response is extremely polarization dependent for superconductors with an anisotropic energy 
gap. Moreover, it was shown that the collective modes which appear in the case of rf— wave 
superconductors are of little importance to the Raman response |^,^ . We can then separate 
the self energy in two parts 

E(q,cj) = S(q = 0,cu) + 5S(q,u;), (5) 

Delaying a discussion of 5S until Section III, we can write down the spectrum of the self 
energy at g = in the pair approximation, e.g., neglecting collective modes as 

Sp,p(q= 0,0;) = {—^ / )tanh(w/4T). (6) 

^ Vcu2 - 4 I A(k) |2 

The subscript g, g denotes the pair susceptibility calculated with vertices g. The real part can 

be obtained via a Kramers-Kronig transformation. Here (...) denotes an average over the 

Fermi surface, Np is the density of states per spin at the Fermi level, G is a Theta function 

and A(k) is the generalized k-dependent energy gap. We see that if the gap is isotropic, 

(A(k) = A), the average around the Fermi surface is frequency independent and thus the 

symmetry of the vertex only determines an overall prefactor of the self energy. Also, since the 

imaginary part of the self-energy has a divergence at the pair threshold energy 2A a phonon 

with a frequency below the threshold should be infinitely sharp (neglecting strong-coupling 

effects). However, if the gap is anisotropic, the vertex and gap couple when averaging over 

the Fermi surface to produce non-trivial changes in the self-energy of phonons of different 

symmetries. Further, if the gap vanishes on the Fermi surface, the presence of the nodes 

can provide decay channels for the phonon leading to a finite linewidth for all non-zero 

frequencies 0]. 

The isotropic {L = 0) density-like terms will be coupled to the long range Coulomb forces 
and thus we must take screening of the vertex into account. Summing R.P.A. diagrams we 
recover the known result at g = 0, 

= ^g,g - (7) 

where 1 denotes the L = contribution of the vertex g |TD|,|TT|. Therefore we see that the 
L = terms are completely screened for q = as a consequence of the long range Coulomb 
interactions and do not contribute to the Raman response. 
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Carrying out the integrations in Eq. (6) using a dj.2_y2 gap, A(k, T) = Aq(T) cos(2(/)) for 
a cylindrical Fermi surface, the spectrum of the phonon self energies can be written down 
in terms of complete elliptical integrals K and E of the first and second kinds, respectively. 
Taking screening into account and defining x = uj/2Aq, we obtain for T = 0, 

Sttx \x[{1 + 2x^)K{1/x)-2{1 + x^)E{1/x)], x > 1, 

(8) 

i.e., the Big channel is not affected by Coulomb screening, while 

with the spectral functions 

r [(7 - 8a;2 + \^x^)K{x) - (7 - \2x'^ + 'i2x^)E{x% x < 1, 

X < (10) 
U^[(32 -28/x2 + ll/x^)K(l/x)-(32-12/x2 + 7/x^)E(l/a;)], x > 1, 



-2y2iV^^^,^ r [(1 + 2x')K{x) - (1 + 4x2)E(x)], x < 1, 

ii„iVH ^^^^^ j (l/x)[(4 - - (4 + x > 1, 



(11) 

and 

Sii(q = 0,a;) = < (12) 

The response functions for finite T are obtained simply by multiplying Eqs. (8) and (10-12) 
by the factor tanh(ci;/4T). The partial screening of the A\g channel by long-range Coulomb 
forces comes from the observation that the square of the energy gap enters into the response 
function in Eq. (6). For the case of d^i_y2 pairing symmetry, the energy gap squared contains 
a term which transforms according to Axg symmetry which leads to a mixing of the L = 
and L = 4 A\g basis functions. This corresponds to partial "transverse screening" of the 
A\g vertex 0. 

The corresponding real parts were obtained via Kramers-Kronig analysis and are plotted 
together with the imaginary parts in Figure 1 for the B\g and screened A\g channels. We 
see that the peak in the imaginary part of the self-energy (which determines the linewidth of 
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the phonon) lies at different frequencies Upeak ~ 2Ao(T) and 1.2Ao(T) for the Big and Aig 
channels, respectively. This is a consequence of the angular averaging which couples the gap 
and e — ph vertex, and leads to constructive (destructive) interference under averaging if the 
vertex and the gap have the same (different) symmetry. Similar behavior for the electronic 
contribution to Raman scattering led to the reasoning that the symmetry which shows the 
highest peak position gives an unique indication of the predominant symmetry of the gap 

0. The symmetry dependence is also manifest in the low frequency behavior, which can be 
written as 

^'k. (^ - 0) = 3NFgiy/A + Oix'), 

^':,Juj-^0)=Npgl^x + Oix'), (13) 

1. e., the spectrum of the self energy rises slower in the Big channel than the Aig channel 



1^ . The power laws are a signature of an energy gap which vanishes on lines on the Fermi 
surface, but the channel dependence of the exponents are unique to a d^2_y2 pair state. 
These channel-dependent power-laws have been observed in the electronic contribution to 
Raman scattering in BSCCO, YBCO, and double and triple layer Thallium cuprates which 
constitute strong evidence for a ci— wave gap of this symmetry as opposed to d^y, d^z or dyz 
symmetry, which also have nodes on lines on the Fermi surface 0. 

The real parts of the self energies (which determines the frequency renormalization) show 
a change of sign near the peak in the imaginary part. While the Big channel shows a mild 
frequency dependence away from the peak maximum and then a rapid change of sign at the 
peak, the Aig channel shows a smooth crossover from negative to positive values, with a 
change of sign that occurs at a frequency which is slightly greater than the peak maximum 
in the imaginary part. Thus a phonon of Aig symmetry which lies at energies below 2Ao(T) 
can become hardened as opposed to softened. We immediately can draw the conclusion 
that phonons of the same frequency will show qualitatively different behavior in different 
channels as a consequence of their symmetry. Therefore, careful attention must be paid to 
symmetry before an analysis of the gap can be made by locating the point where phonon 
softening or hardening occurs. 

III. TEMPERATURE DEPENDENCE 

We now investigate the temperature dependence of the phonon lineshape. The q = 
spectral function, Eq. (6), vanishes at due to the lack of particle- hole continuum for 
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pair creation. This term thus always predicts phonon broadening compared to the normal 
state below Tc for a gap with nodes. However, the term responsible for the normal metal 
self energy (due to, eg., finite momentum transfer or anharmonic decay) will be affected by 
superconductivity due to the reorganization of the density of states as the gap opens up. 
In order to recover the normal metal lineshape at Tc, one must use finite q (or impurity 



scattering |JT0|,|13[) to generate the additional term 5Z1 which does not vanish at Tc. We now 
generalize the result to finite q for anisotropic gaps. For finite g, the spectrum of 5S at finite 
temperatures is given by 5S"(g,T, cj) = Qivpq — ^^)(| g]^ P -^(k, cj)) with 

2vFq ^|A(k)| ^ Je2_ I ^(k) \^J{E + uy- I A(k) |2 



+eiE- I A(k) I -umE-u)-fm EiE-.y\Aik)\^ ^^^^ 

Ie^- I A(k) \''J{E-ujf- I A(k) 



where / is a Fermi function. The Theta function B(fi?g — c<j) restricts the frequency shift due 
to phase space consideration, reflecting that the region of the particle-hole continuum van- 
ishes for small wavenumbers as a consequence of momentum conservation. Since fj^g <^ A 
in the cuprate materials and also in A- 15 materials, this term will only contribute to the self 
energy for phonons of very small energy. However, it has been shown for s— wave supercon- 
ductors [jlOl that the incorporation of impurity scattering removes the phase-space restriction 
due to the lifting of momentum conservation and (5E contributes for all frequencies. While 
incorporating impurity scattering remains beyond the scope of the present treatment, we re- 
mark that it is expected that a similar consideration for the case of (i— wave superconductors 
would also lead to the contribution of (5S for all frequencies. This remains to be explored 

0- 

In the limit of small frequencies (cj << T), we obtain the simple result 

5S"(g, uj«T)= j^^j' )Q{v,q - u). (15) 

vpq el'^wl/^ + 1 

Similarly, at Tc, Eq. (14) recovers 5'Z"{q,uj « T^) = Q{vpq — uj)uj^^^^^{\ P) and thus the 



Vpq \l .^k 

ratio of the low frequency response in the superconducting state to that of a normal metal 



at Tc is given by 

5S"(g,cu « T) 





V(|A(k)|)) 







(16) 



6T."{q,uo « Tc) 

This shows how the redistribution of the density of states below Tc to higher energies as the 
gap opens up leads to a reduction of the decay channels available to particle-hole creation and 
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a net decrease in the phonon linewidth. In isotropic superconductors, the Fermi function 
can be pulled out of the average and the resulting expression is independent of phonon 
symmetry. However, once the gap is anisotropic, there exists coupling between the vertex 
and the the gap which leads to a symmetry dependent result. 

Using a weak coupling expression for the temperature dependence of the energy gap 
(2Ao/fcsTc = 5.1252), we numerically evaluate Eq. (16) while taking screening into account. 
The results are plotted in Fig. 2 as a function of T/Tc for a dj.2_y2 energy gap compared 
to a BCS isotropic gap. The low temperature behavior is given by a power-law in T for all 
channels for the (i— wave case while the ubiquitous exponential dependence in T is seen for 
all channels in the s— wave case. The power-law behavior for the ci— wave case is channel 
dependent, with exponents identical to those of Eq. (13), in the sense that u can be replaced 
by T. What is remarkable is that the fall off of the Fermi function at low temperatures is 
quite slow in those channels which are orthogonal to the symmetry of the gap, with the 
notable example of the Aig channel, which shows a residual broadening at T/Tc = 0.3 of 
roughly 20 percent of that of the normal state. This was argued in the case of electronic 
Raman scattering to be further evidence for an energy gap in the cuprate materials which 
has predominantly Big character, due to the observation that a gap opens up quickly in the 
Big channel compared to Aig and others which have been probed via Raman scattering [0 . 

For higher frequencies > T, we have evaluated Eq. 14 directly. The results are 
quite similar to those of Fig. 2 for all frequencies uj up to roughly 4A, but then at higher 
frequencies all channels eventually display a linear T dependence (i.e., the behavior of the 
normal state) for energy scales much greater than the gap energy. 



IV. ENTIRE SPECTRAL FUNCTION AND ROLE OF SATELLITES 



In this section we consider the entire phonon spectral function Eq. (1), paying particular 
attention to the role of satellites which arise due to e — ph coupling. The role of satellites 
have not been explored in anisotropic superconductors. As is well known for BCS supercon- 
ductors, satellites appear in the phonon spectral function for all frequencies of the optical 
phonon, but have the greatest residue for phonons near twice the gap edge. In the BCS case, 
impurities wipe out the satellite peak [jlO[, explaining why they have yet to be definitively 
observed in conventional A-15 superconductors. In the absence of impurities however vastly 
different lineshapes can be obtained due to the interference of the satellites. 

Using a gap of dx2_y2 symmetry and working specifically at g = 0, Eq. (6), we find 
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that the satelhtes are present in anisotropic superconductors as well due to the fact that 
the real part of the denominator of Eq. (1) has two zeroes for any value of c^o - one at the 
renormalized phonon frequency and the other at the satellite position. The satellite becomes 
more pronounced the closer the optical phonon is to the peak position of the self energy 
spectrum as in the BCS case and interferes with the phonon. Therefore, for the case of a 
phonon located below the spectral maximum, where the satellite peak is observable only 
at T=0 for large e — ph coupling, as the gap decreases on approaching Tc the satellite will 
be made to pass through the phonon position and will be subsequently distorted. This is 
shown in Fig. 3 for a Big phonon (cuo/Ao(r = 0) = 1.0, c/|j^^A^f/Ao(T = 0) = 0.1) for the 
temperatures indicated. The phonon lineshape is drastically affected by the satellite which 
takes spectral weight away from the phonon when the peak of the spectral function is close 
to the phonon position. The phonon linewidth grows as the peak of the spectrum moves 
up in energy with decreasing T and is hardened. The linewidth and frequency shift reaches 
a maximum when the peak and phonon position coincide and then the linewidth decreases 
and the phonon softens as T — Tc. 



V. CONCLUSIONS AND COMPARISON WITH EXPERIMENT 

In this Section we combine the previous results and examine the phonon linewidth as 
a function of temperature for the case of two phonons which lie at approximately 285, 340 
cm~^ for the Big channel and 464,500 cm^^ for the Aig channel in BSCCO and YBCO, 
respectively. Using our previous fits to the electronic Raman scattering in BSCCO, we 
obtained a value of the energy gap at T = to be Ao(T = 0) = 287 cm~^ [0. Therefore the 
normalized optical phonon frequency is given by u;o/Ao(T = 0) ~ 0.99, 1.62 for the Big, Aig 
phonon, respectively in BSCCO, while for YBCO the ratio is expected to be slightly higher. 
We can immediately make the following statement. Since the interference effects of the 
phonon with the satellite peak can only occur for a phonon which is located at T = 
below the peak in the imaginary part of the self energy, there should be no interference 
effects on the Aig phonon since it lies above the peak in the spectrum at T = 0. Thus its 
renormalization should be a monotonic function of temperature. However that is not the 
case for the Big phonon. Anomalous behavior of the Big renormalization as seen in Fig. 
3 arises due to the interference between the phonon and the rapid rise of the self energy 
near 2Ao(T), which passes through the phonon frequency at T/Tc ~ 0.9. Another remark 
is in order. In order to make an accurate fit to the data, the magnitude of the coupling 
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constant needs to be addressed. As we have seen in Section IV, it controls the strength of 
the satelhte and its subsequent effect on the phonon hneshape. Little is known about the 
coupling constant ||^ and thus we can only make general statements on the behavior of the 
phonons. The magnitude of the effect cannot be predicted. 

Inspecting Fig. 1, the q = part of the self energy, at each phonon frequency, we see 
that the Big phonon is broadened and softened at T = compared to while the Aig 
phonon is broadened but lies right at the point where the real part is changing sign. This 
term most accurately describes what is seen in the phonons in YBCO. A rapid rise of the 
Big phonon linewidth below has been seen |]6|,[T5|, reflecting the interference of the peak 
in the self energy and the phonon, (see Eq. (6) and Fig. 3). The Big phonon additionally 
is mostly softened to lower frequencies. For the Aig phonon, the frequency shift is seen only 
to slightly higher frequencies. This is due to the real part of the self energy crossing zero at 
the phonon position. However, the Aig phonon narrows just below Tc, reaches a minimum 
width and then broadens at lower temperatures. This could be a result of a competition 
between the effects of S and 5S, but without information concerning the magnitude of the 
coupling constant, finite momentum transfer or impurity scattering this remains an open 
question. 

In addition, we have the contribution arising from non-zero g, Eq. (14), which indicates 
that this contribution to the self energy is reduced compared to its value at T^. and decreases 
as and T for low temperatures in the Big and Aig channels, respectively (see Fig. 2). 
This term most accurately describes the experiments in BSCCO. The linewidth of both the 
Big and Aig phonons decrease monotonically with temperature, which points to the lack of 
contribution coming from satellite effects. The linewidth decrease in the Big channel can 
be fit with a dependence while a term linear in T can be fit to the Aig phonon which is 
the same dependence as in the normal state. Both behaviors are consistent with Eq. (16). 
This is to be compared with the predictions of an isotropic s— wave theory lines, which are 
identical for each channel (~ e~^/"^). The theory for a gap of Big character shows a marked 
symmetry dependence resulting from the interplay of gap and vertex symmetry. While it 
is arguable whether an exponential temperature dependence can also fit the Big data [^, 
the lack of change of the exponent of the Aig phonon's linewidth is a direct consequence of 
the energy gap anisotropy. Therefore it will appear that the Aig phonon will be unaffected 
by superconductivity. Since the Big phonon has the same symmetry of a d^i-yi gap, it's 
linewidth will show the greatest change due to the onset of superconductivity. 

In conclusion, within the accuracy of the experiments, we have seen that the changes 
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in the phonon lineshape as a function of symmetry can be explained with a choice of the 
gap which has (at least predominantly) Big character, supporting recent comparisons made 
on electronic Raman scattering on BSCCO. However, without knowledge of the magnitude 
of the coupling constant, the importance of finite q and satellite effects remains an open 
question. Of course other choices of gaps which have a small but finite minimum value, eg. 
anisotropic s-wave or s + id would give similar results to the dx2-y2 choice for the gap. Both 
the electronic and phonon contributions to Raman scattering below Tc can be explained 
by simply invoking the symmetry of the vertex which couples to the symmetry of the gap. 
More detailed experiments would be extremely useful to pin down the magnitude of the 
e — ph vertex and subsequently the role of the satellites, and the role of impurities and the 
mechanism and magnitude of the coupling remains to be explored 
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FIGURES 

FIG. 1. Real and imaginary parts of the phonon self energy for the Big and Aig channels for 
a cylindrical Fermi surface. Magnitude of the vertices are set equal to one. 

FIG. 2. Temperature dependence of the a; — imaginary part of the self energy in a rf^^-j/Z 
and isotropic BCS superconductor compared to the normal state. 

FIG. 3. Phonon spectral function for a Big,{uo = Ao{T = 0),gl^^NF//^o{T = 0) = 0.1), 
phonon in a superconductor with dj.2_y2 pairing symmetry for various values of T/Tc as indicated 
in upper part of Figure. 
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